Do non-relativistic neutrinos constitute the dark matter? 
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PACS 95 . 35 . +d - Dark matter 

PACS 98 . 65 . -r - Galaxy groups, clusters, and superclusters 

PACS 14. 60. St - Non-standard-model neutrinos, right-handed neutrinos, etc. 

Abstract. - The dark matter of the Abell 1689 galaxy cluster is modeled by thermal, non- 
relativistic gravitating fermions and its galaxies and X-ray gas by isothermal distributions. A fit 
yields a mass of h^g^ {12/ g^^^LUb (30) eV. A dark matter fraction = /i^^f ''^0.1893 (39) occurs 
for g = 12 degrees of freedom, i. e., for 3 families of left plus right handed neutrinos with masses 
« 2^^'^Gy^m^. Given a temperature of 0.045 K and a de Broglie length of 0.20 mm, they establish 
a quantum structure of several million light years across, the largest known in the Universe. The 
virial a-particle temperature of 9.9 ± 1.1 keV/fcs coincides with the average one of X-rays. The 
results are compatible with neutrino genesis, nucleosynthesis and free streaming. The neutrinos 
condense on the cluster at redshift z ^ 28, thereby causing reionization of the intracluster gas 
without assistance of heavy stars. The baryons are poor tracers of the dark matter density. 



Introduction. — Dark matter is postulated by Oort 
to explain the motion and density of stars perpendicular 
to the Galactic plane [1]. Zwicky points out that galaxy 
clusters must contain dark matter, [2] while Rubin demon- 
strates that galaxies require dark matter in order to ex- 
plain the rotation curves of stars and hydrogen clouds [3] . 
Nowadays, gravitational lensing observation is standard- 
ized, and dark matter filaments on the scale of clusters of 
galaxies between empty voids can be inferred [4] . 

The main dark matter candidates are Massive Astro- 
physical Compact Halo Objects (MACHOs), Weakly In- 
teracting Massive Particles (WIMPs), elementary parti- 
cles, and further, e. g., axions. Of the total mass of the 
universe, Hb — « 4.23% consists of baryons, 

of which a minor part is luminous, a part is located in 
gas clouds, and a part makes up the galactic dark matter 
as MACHOs [5-7]. However, the remaining mass of the 
universe is non-baryonic, dark energy and dark matter. 
The dark matter, with cosmic fraction ilu = 21.4 ± 2.7% 
according to WMAP5 [8], will be the focus of the present 
work, where we assume it to consist of fermionic WIMPs. 

Thermal fermion model. — We consider non- 
interacting dark fermions (x) with mass m and g degrees of 
freedom, subject to a spherically symmetric gravitational 
potential U (r) and in equilibrium at temperature T. The 
mass density reads 



J (27r?i)3 exp[(p2/2m + mC/(r) - fiJ/keT] + 1 ' ^ ^ 

where p is the momentum and fi = aksT the chemical 
potential. The spherically symmetric potential normalized 
at U{0) = reads f/(r) = Gj d^r' p{r'){l/r' - l/|r - r'|). 
It satisfies the Poisson equation U" + 2U' /r = AnGp. In 
dimensionless variables, x = r/R^, (p = mU/ksT, one has 



Px[r) = -^Li3/2 



1 — 



(2) 



with the polylogarithm Li-y(z) =J2T=i^''/^'^ ^'^^ \^\ ^ ^ 
and analytically continued elsewhere. The thermal wave- 
length At and a characteristic scale i?, read, respectively, 
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Accounting also for the galaxies (G) and the hot X-ray gas 
{g), see below, the Poisson equation will take the form 



-0' = -Li3/2 (-e"-^)+e"^ 
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For a — CX3 it reduces to a two component isother- 
mal model, that describes galaxy rotation curves well; for 



-oo it becomes the fermionic isothermal model. 
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The total mass inside a sphere of radius r — R^,x reads 
ksTR 
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Abell 1689. — This is the best studied galaxy cluster, 
known for its large lensing arcs, notably one at the Ein- 
stein radius rE — 50". It is well relaxed and spherically 
symmetric, with an intruding subcluster in the North- 
East, which does not affect the South-West hemisphere. 
It offers a test for the above. From gravitational lensing 
the azimuthally averaged mass profile is deduced, [9] 
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The average redshift is z = 0.183. Observations are 
presented up to — h^^ Mpc, with Hubble constant 
Hq = lOO/i km/s Mpc. Relatively small fluctuations occur 
in the contrast function AE(r) of the averaged E within 
r versus outside r. The first piece reads 



= A / drx2^ri S(ri) = -l-M^oir). (7) 
7rr^ Jq 7rr^ 

The average between r and is directly related to this, 
S(r^r„) = [M2D{r^) - M2D{r)]/[nirl - r^)], so 



AE(r) = E(r)-S(r^r„0 = 
From (7) and (6) one derives 



S(r) - E(r^) 
1 - r2/r^ 



(8) 



E(r) ^ — / dr'r'p{r') 



(9) 



After defining an amplitude A = hP /2g'^G^m^ R!1, elimi- 
nating T via (3), T = Trh'^/2fG^kBm'^Rt, and using (4), 
also S can be expressed in terms of 



ds0'(xcoshs). (10) 



To proceed, we consider the baryonic matter, galax- 
ies (G) and, mostly, a hot X-ray gas (g). In hydrostatic 
equihbrium, p'J pi = —GM{r)/r'^, {i = x,G,g), both the 
classical galaxies and the low density gas have a Boltz- 
mann distribution, pQ ~ exp(— J7/cr^) with ay the line- 
of-sight velocity dispersion and pg ^ eiq){—l3grngU) with 
Wig = O.GOQrriN the average mass in the gas.^ Hydrostatic 
equilibrium then imposes f3g = l/ksTg where Tg is the gas 
temperature. With I3g = ksT Imaf, and (3g = WigT/mTg, 



^For typical 0.3 solar metallicity one has njj = lOriHci so pg 
2.3niikgTg, pg = lAniirrif^ and mg/mpf = 1.4/2.3. 



Fig. 1: The mass contrast AE as function of radius r. Large 
data points from Ref. [9], small ones (at radii 5 2"^^h^g kpc 
with n = 0, ..,12) constructed from Fig. 6 of Ref. [10]. Full 
line: The theoretical profile; it softens below 7 kpc. 



the densities may be written as pi = gmX^^ exp(Q!i — Picj)), 
{i = G,g). This leads to the last two terms of Eq. (4). 

We can now make a fit of Eqs. (8) and (10) to the 19 
data points of [9], combined with 13 points constructed 
from recent core values for M2d{i') [10] and S(rm) from 
our model. As seen in Fig. 1, the relative errors increase 
strongly with r, which reduces the effective number of 
points. The errors become more equal if we consider the 



As explained later, we take (3g = 0.153 and 
ag — 2.36 at Pg — 1- There is a minimum — 13.617 at 
Pg — 0.80. This is close to the virialized value /3g = Ij so 
we stick to that with its = 13.645. With h = 0.70 /170, 
the correlation matrix for the upper errors yields 



^ = 59.4 ±9.6 /l7oM0pc-^ a = 38.4 ±3.1, 

i?* = 297± 10/ifo^kpc, ac = 8.26 ±0.32. (11) 



We present its fit in Figure 1. The WIMP mass reads 



2l/8gl/4 ^3/8^1/8^5/8- 



(12) 



Since ARl = 136 ± 25 h^^M^G^C', m has a 2% error. 
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1.455 ±0.030 eV. 



(13) 



With T-yo = 2.725 K, the global fcrmion density is 
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While g is the number of states that can be filled in the 
cluster formation process, g is the filling factor in the dark 
matter genesis. The global mass fraction thus reads 
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Neutrino mass and neutrino dark matter 
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= /lyo^'^O.lSgSi 0.0039. (15) 
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The gravitino of supersymmetry (s=3/2, g — 8) can 
explain the Abell data, but it decouples early, in the pres- 
ence of ~ 100 relativistic degrees of freedom, so that 
g ~ 0.4 leads to a small ft^ ~ 0.8%. The same holds for 
other early decouplers [11]. Bosons, like the axion, can 
not fit the data because of the tilt in Fig. 1 at r < 200 
kpc; axionic Bose-Einstein condensation can only exist up 
to the small scale yOvrR*, which is of no help. 

Neutrinos. — They can occupy in the cluster forma- 
tion process all g = 12 left and righthanded states, which 

1/2 

gives m ~ 1.455(30) /ijq eV. Neutrinos oscillate, [12] 
Am?2 = S.Oto glO-^eV^ Am^g = 1.9 - 3.0 IQ-^eV^ so 
the mass eigenvalues differ. It is natural to suppose that 
all virial speeds are equal, so Ti — Trrii/m. If also their 
chemical potentials behave as fii — ksTia, the above ap- 
proach still applies with m = (mi -I- 7712 + m3)/3. The 
degeneracy parameter will change negligibly, 
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The oscillations bring a shift smaller than our error bars. 
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Sterile masses are expected to weigh keV's or more, see 
[13] for a review, but a (near) equality between left and 
right handed masses and their abundances is needed to 

r-i i,^3/2n 19 v~>12 // 1 v~>12 4M/4 

maximize il^ w /lyp ^ l^i=i nim.^/ [j^ l^i^i nim^y^. 

The de Broglie length in the cluster \tu — 
2gn-^Gm^Rl = 0.20 mm is visible to the human eye; the 
Compton length is 0.136 /j,. The A1689 neutrino temper- 
ature = TTGm^AR^:/kB — 0.0447 K is low and makes 
them strongly non-relativistic, with local dispersion 
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Fig. 2: Blue line: The neutrino number per cubic thermal 
length and per degree of freedom as function r. The point 
(505 /lyQ^kpc, 1) separates the region A'^^ > 1 which exhibits 
strong quantum effects from the classical region A^^ <C 1. 
Green line: The galaxy mass density, at the same scale, is 
concentrated in the center. Red line: The gas mass density at 
the same scale. 

0.34 GeV/cc, while ns(0) = 5.1 10^/cc. The quantum pa- 
rameter Nyijr) = — Li3/2(— e"^"^) is plotted in Fig. 2. Its 
maximum is N^[0) — 180, so quantum statistics is indis- 
pensable. We may define the quantum-to-classical tran- 
sition by N^{rqc) = 1. This gives r^c = 505/ifo^kpc or 
diameter 3.3 10^ lyr, a giant size for quantum behavior. 

As seen in Fig. 2, the baryons are poor tracers of the 
dark matter density, even they do trace the enclosed mass. 

Neutrino free streaming (fs) in expanding space, H{t) 
= a/ a, is described by a coUisionless Boltzmann equation. 



dtfu H \VH 



^Gm^M{r, t) 



dpU = 0. (19) 



Below the Compton temperature = 16, 850 K the dis- 
tribution fl^ ~ \l(^gPc/kBT^ _|_ 2) is long maintained with 
shding Ty = (4/ll)^/^T-y, until the last, Newtonian term 
sets in. In the fs regime the scaling p ~ ksT^/c implies 



Pfs 



ypy^^ 2.567^. 



(20) 
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which at large r equals the galaxy velocity dispersion = 
y/kBT/m — 488 ± 60 km/s. The latter agrees reasonably 
with estimated speeds in [10] and with the 295 ±40 km/s 
of the singular isothermal sphere that fits the mean galaxy 
distribution [14]. Indeed, between 5 and 80 kpc pa looks 
somewhat like a singular isothermal distribution. 

Neutrinos are abundant, see (14), but their speed is too 
low to leave traces such as Cherenkov radiation. Con- 
densed in clusters, their local density is large. One has 
n^(0) = -gX~^U3/2{-e°') = 2.3 lOVcc and p^(0)c2 = 



The neutrinos condense (pc) on the cluster when the typ- 
ical Newton force Gm^M{R^)/Rl = 6.5610 



-10, 



becomes comparable to the free streaming force, ptsH — 
2.82 lO-^^mt,{z + l)y/nA + ^m{z + Ifms-'^, which hap- 
pens for flM ^l-^A^ 0.25 at z^c = 28.3, T^'^ = 77.2 K 
and age 120 Myr. Consistency occurs by the near match 
of kinetic energies, ^■m^[a';^{R^)]'^ = 0.64 x {p'^^)'^ /2m^. 
The free streaming stops then and the temperature sinks 
only moderately, from the cross over value T^""" = 
|m^(crf)VfcB = 0.13 K to present one, = 0.0447 
K. Literature often mentions 1.95 K as present neutrino 
temperature; this derives from relativistic free streaming, 
but since they are non-relativistic, the kinetic tempera- 
ture would even be lower, T^'" = 2E^^.^/5kB = 9.6 lO^-^K. 
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Our is higher because the condensation basically stops 
the cooling. The neutrino gas then clusters with a speed 
of sound Wg'^ — yjbp/ip = 2620 km/s at the Jeans 
scale Lji, ~ v""^ j ^fG(^ = 1.57 Mpc. The Jeans mass 
Mj^ = ttpI'^L'^j^/Q = 1.3 lO^'^AifQ estimates the total mass 
(5) at Tm = lMpc//i, Mtot(r„) = 6.91O"M0. Only 1.1 
% of this is in galaxies, and 1.5% in gas: the A1689 cluster 
is baryon poor, a property noticed before [9, 16].^ 

Virialization, gas profile and X-ray emission. — 

Partial thermalization takes place since during the conden- 
sation the individual objects, neutrinos, H, He and other 
gas atoms, and galaxies move in a time dependent grav- 
itational potential, which changes each ones energy at a 
rate proportional to its mass. If the phase space occupa- 
tion becomes uniform, a Fermi-Dirac distribution emerges 
[18] and our approach applies; else it yields an estimate. 

This thermalization heats and reionizes the loosely 
bound parts of the intracluster gas. The virial proton tem- 
perature of the A1689 cluster, = nipT^/m^ = 2.48 
keV, has the order of magnitude of the average X-ray 
temperature 10.5 ± 0.1 keV of the cluster and the one of 
the South Western hemisphere, 11.1 ± 0.6 keV [15]. But 
virial equilibrium of the a-particles, at temperature = 
maT^/mv = 9.9 ± 1.1 keV, explains these observations. 
For this reason, we have already taken f3g = rn.g/ma- 

As expected theoretically, this gives a good match for 
the mass profile of the gas, which is deduced from the 
X-ray profile [16]. Fitting their first data point Mg = 
7.79 h~^MQ at 30ft.~^kpc to our results sets ag = 2.36 at 
h = 0.7, as used above. Then the next 6 data points, up to 
200 h~^kpc, match within the symbol size with our theory. 
From there on, it underestimates the data, by a factor 4.5 
at the last point at 716ft~^kpc. Reasons for this may be 
a decay of the temperature beyond 200/i~^kpc [16] and a 
contamination of the data from the more relaxed South 
West with those from the less relaxed North East. 

The bolometric (total) energy emission 1.66 ± 0.64 10"^^ 
W [16] must be supplied by a slow contraction of the 
cluster. Thus the radiation serves to maintain the virial 
equilibrium state as it does for stars [19]. This will double 
the A1689 gravitational energy in some 700 Gyr. 

Extended standard model and dark matter frac- 
tion. — Righthanded neutrinos have no hypercharge, so 

they justly do not enter the Z decay. Their mass is de- 
scribed as for quarks. The Yukawa eigenvalue Y,^-^ = 
can be expressed in the electron coupling 
Ye = 2.935 10-^ as Y^^ = 0.966(20) hl(,^Y^. This exhibits 
some order in the lepton masses; in case hro = 1.062(42) 
orh = 0.744(30) it impUes 



m 



= YeTUe = 2^/'^G^J^ml = 1.4998 eV. 



(21) 



^The missing baryons may be located in intcrcluster clouds, such 
as the observed 0.91 kcV gas of mass ~ 10" M0 in a bridge between 
the 4.43 keV A222 cluster and the 5.31 keV A223 cluster [17]. 



Active neutrinos (lefthanded neutrinos, righthanded an- 
tineutrinos) have .9 = 6 degrees of freedom. Eq. (15) leads 
to a cosmic density 0.0952 ±0.0019 h^Q^^, clearly exceed- 
ing the 0.028 /ifcf of WMAP5 (0.013 /if^f when combined 
with baryon acoustic oscillations and supernovas) [8]. 

The situation can be even more interesting, since the oc- 
cupation of the (mostly) righthanded states (sterile neu- 
trinos) can have become sizeable if there is also a Ma- 
jorana mass matrix. The latter couples the neutrino to 
its charge conjugate, rather than to the antineutrino, see 
e. g. [20]. This allows neutrinolcss double /3-decay, 
where the two neutrinos emitted in the /3 decays anni- 
hilate each other, and the electrons leave with opposite 
momenta. For simplicity we consider the case of 6 ster- 
ile states, bringing the total number of neutrino states 
at g = 12, the case discussed so far. In order to keep 
nearly equal masses, we need small Majorana terms Mj 
{i = 1,2,3 denotes the families), not the large ones of 
the see-saw mechanism. A one-family version of the prob- 
lem shows that a (nearly) thermal occupation of sterile 
modes is possible above the decoupling temperature 3.5 
MeV, provided M, > ZlQ-^eV'^/m^ « 2 lO'^eV [21]. 
Experimental searches have determined the upper bound 
« rrifff} < 0.2 - 0.7 eV [22]. This filling implies that 
g = 12 can be reached, leading to a dark matter fraction 



fl^ < 0.1904 ±0.0038/1, 
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The case g g > 12 is also possible. For 5 ~ 33 it has 
enough matter to reach to ~ 1.0 /Hq without dark energy. 

Nucleosynthesis. — Our additional relativistic mat- 
ter can be coded in the enhancement factor S [23], p' = 
S^p; after e"*" — e~ freeze out it reads for g = 12 



16 + 7£f(4/ll)4/3 

16 + 42(4/11)4/3 
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= 1.1854. 



(23) 



This enhances expansion, H' = SH, leaving less time for 
neutron decay and resulting in too much ^He. It can be 
balanced by a neutrino asymmetry due to a dimension- 
less chemical potential ^, an effect which induces more n 
decays via n + Vg ^ P + e, 



Le = 



12C(3) \T. 



T 



With 7710 = IO^^hbI 



1210^/170 J the cosmic mi- 



crowave background (CMB) value reads T]y^^ — 5.60 ± 
0.15 [8]. The "^He value 0.881^;^^ does not fit to it, which 
motivated to put conservative error bars [23]. But the 
effects of extra matter and asymmetry on He are large, 
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So we may fix ^ by matching to i]^^^ . This yields ^ = 
0.162, Le = 0.040. Other authors report a ^Hc value closer 
to the one of CMB [24] , but this does not modify ^ much. 
For Li, TjLi = 6.051q ^[2: 9^ similar approach brings [23] 

mo = rjLi + 3(5 - 1) + = 6.88t!];}3, (26) 
while for deuterium, rin = 5.92+q;33, it implies 

mo=»?D + 6(5-l)-^^ = 6.83l°:|°. (27) 

The freedom in ^ appears to solve the *He discrepancy. 
Since the CMB value still has to be rcdcrivcd for the prior 
of neutrino dark matter, the final result may end up near 
rjio = 6.88 ± 0.15, riB^o = 5-69 ± 0.12%, to be compared 
with Ob/?,?(, = 4.63 ±0.12% from WMAP5. Together with 
(22) it would lead to a total matter fraction I^m = i^B + 
0^ w 24.7 ± 0.5%, while WMAP5 reports 25.8 ± 3.0%. 

Conclusion. — On the basis of three assumptions, 
Newton's law, quantum statistics and virialization, we de- 
rive the profile of quantum particles (WIMPs), galaxies 
and intraclustcr gas. Because of the virialization — equal 
velocity dispersions for WIMPs, galaxies and a-particles, 
or, more precisely, each ones temperature proportional to 
its mass the WIMPs have a polylogarithmic profile, 
while the galaxies and the gas have isothermal profiles. 

A fit to total (lensing) mass observations of the cluster 
Abell 1689 is possible only for fermionic WIMPs with eV 
mass. The error of this method is small, 2% for the present 
data set. Although we have not shown the validity of our 
virial equilibrium assumption, the explanation of the X- 
ray temperature of the hot gas, Tg ~ lOkeV = 116 lO^K, 
as the virial temperature of a-particles is striking. Due 
to collisions that temperature is shared by the electrons, 
protons and ions. The predicted mass profiles for galaxies 
and gas are also consistent with observations. 

This suggests that our approach cannot be far off, so the 
WIMP mass is a few eV and dark matter is hot. Early 
decouplers that have been in equilibrium would yield a 
small cosmic dark matter fraction [11], so if they would 
set the dark matter of the A1689 cluster, there should also 
be other dark matter, the major part, but absent in this 
cluster, which is unlikely. Therefore early decouplers such 
as the gravitino are ruled out as dark matter candidates. 
The thermal axion is ruled out because it is a boson. 

The case of g = 12 degrees of freedom performs well, 
pointing at three families of left and right handed fermions 
and antifermions. The obvious candidate is the massive 
neutrino, because when condensed in the cluster the left 
and righthanded states are equally available. The mass is 

1 /2 

then — 1.445 h-jQ eV with a 2% margin and smaller 
variations between the species due to neutrino oscillations. 
There is the striking connection m,^ ^ The 
dark matter fraction of active neutrinos is then 9.5%, 



showing that the cold dark matter analysis, that allows 
only 1.3% at best [8], must definitively be erroneous. 

The scenario in which dark matter is, say, half due to 
neutrinos and half due to cold dark matter (CDM) par- 
ticles was found viable in connection with violent relax- 
ation [25]. But it does not fit the A1689 cluster. Indeed, 
heavy particles have a Boltzmann distribution. Being col- 
lisionless and relaxed, they are accounted for already by 
the isothermal galaxy term in Eq. (4) with /3g = 1. So at 
best they present 1-3% of the A1689 mass, too little for 
this 50-50 assumption, so it would again imply the unlikely 
conclusion that this cluster is not representative. 

The Tremaine-Gunn argument [26] of no increase of 
the maximal phase space density (except for a factor 2) is 
automatically satisfied by the Fermi-Dirac distribution. ^ 

To describe the dark matter profile of a relaxed galaxy 
cluster is a clean problem that involves almost no cosmol- 
ogy, so its daring predictions pose a firm confrontation 
to conclusions based on more intricate cosmological theo- 
ries, such as the cold dark matter model with cosmologi- 
cal constant (ACDM model, concordance model). Indeed, 
our findings are in sharp contradiction with present cos- 
mological understanding, where neutrinos are believed to 
be ruled out as major dark matter source [13]. Studies like 
WMAP5 arrive at bounds of the type m,^^ + rriv^ + m^^ < 
0.5 eV. They start from the CDM paradigm, or from a 
mixture of CDM and neutrinos, the reason for this being 
indirect, namely that without CDM the cosmic microwave 
background peaks have found no explanation. But the 
CDM particle has not been detected, so other paradigms, 
such as neutrino dark matter, cannot be dismissed at fore- 
hand. The CDM assumption has already questioned [5,27] 
and it is also concluded that WIMP dark matter has an 
eVmass [28]. 

The common assumption that light (baryons, Lyman— a 
forest aspects) or intracluster gas trace the local dark mat- 
ter density appears to be invalid, even in the absence of a 
temperature gradient, for radii at least up to 1.5 Mpc, see 
fig. 2. Nevertheless, they do trace the enclosed total mass, 
since this is imposed by hydrostatic equilibrium. The 
galaxies (G) behave differently from the neutrinos even for 
= 1) because the first ones are classical while the sec- 
ond ones are degenerate fermions; though non-degenerate, 
the gas (g) behaves differently from both of them, because 
it is an ionic mixture with electrons, implying a parameter 
< pG- This non sequitur nullifies many conclusions in 
literature, notably that sterile neutrinos should have keV 
mass at least and connections between the Lyman- a forest 
and local dark matter densities [10,13,29]. 

^Let us answer a criticism often met in literature. Galaxies and 
dwarf galaxies may have their baryonic dark matter in the form of 
MACHOs, H-He planets of earth mass [5], thousands of which have 
been observed [6, 7] . In that scenario the (dwarf) galaxy Tremaine- 
Gunn bound involves the proton mass and is satisfied a million times. 

''CMB peaks may arise without CDM seeds, namely from viscous 
instabilities in the baryonic plasma before and at decoupling [5,7]. 
Free streaming WIMPs would then not have time enough to wash 
out the newly created baryonic structures. 
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As mentioned, active neutrinos alone would bring about 
half of the expected dark matter. In the early Universe the 
sterile ones can be created too, at temperatures between 
200 and 3.5 MeV, provided neutrinoless double /3-decay is 
possible, a process which violates the lepton number. Neu- 
trinos then are Majorana particles. The related Majorana 
mass should exceed rn/3/3 > 10~^eV or so, which gracefully 
respects the experimental upper bound m/3/3 < 0.2 — 0.7 

eV [22]. Then < 0.19/?,^^^ 

Both the neutrinos and the antineutrinos fall down- 
wardly in a gravitational field, as usual. 

The neutrinos stream freely, cool in expanding space 
and condense on clusters at z ~ 28, forming their dark 
matter. They simultaneously heat the intracluster gas, 
which ionizes on its way to the 10 keV virial equilibrium 
temperature. This makes it plausible that it is the neu- 
trino condensation on clusters that, as a cosmic virial im- 
print, reionizes all the loosely bound gas, without any need 
for heavy stars, the currently assumed cause [30]. 

The central parts of clusters constitute quantum particle 
structures of several million light years across. We expect 
that the Universe does not contain larger ones, though we 
make a reservation for the unknown cause of dark energy. 

As usual for clusters, we did not have to invoke a modi- 
fied Newton dynamics like MOND [31]. The popular NFW 
mass profile [32] plays no role, since it deals with heavy, 
non-degenerate dark matter, rather than with degenerate 
fermions. Fig. 2 shows that the NFW cusp is actually pro- 
duced by galaxies, classical objects indeed, but not dark. 

The cold dark matter paradigm has to be abandoned, a 
conclusion that was already inferred from observed corre- 
lations in galaxy structures [27]. Hot neutrino dark mat- 
ter has to be reconsidered. A gravitational hydrodynamics 
theory of top - down large scale structure formation has 
been proposed more than a decade ago [5] . 

Our theoretical description of virialized galaxy clusters 
can be used to process observation data. It can be general- 
ized to hydrostatic equilibrium with temperature profiles. 
Precise observations of the profiles of the weak lensing, the 
galaxy velocity dispersion and the X-rays in the relaxed 
South West region of the Abell 1689 galaxy cluster are 
welcome. The prediction for the neutrino mass already 
has a small error and offers, once it is known, a new way 
to determine the Hubble constant. Simulation of the con- 
densation out of free steaming will test our findings. 

As for any physical prediction, the ultimate proof is 
a direct observation, in our case of the neutrino mass. 
The Mainz- Troitsk tritium /3-decay experiment records the 
maximal electron energy in the reaction ^H ''Hc^ + 
e~ + Vei and sets a bound on the electron-antineutrino 
mass, m^^ < 2 eV at the 95% confidence level [33]. This 
leaves our case viable. The KArlsruhe TRItium Neutrino 
experiment (KATRIN), scheduled for between 2012 and 
2015, will test our prediction by searching for a mass down 
to 0.2 eV [33], so that our ~ 1.5 eV regime should be 
relatively easy. A next challenge is to settle the Majorana 



mass involved in neutrinoless double /3-decay. 

The present status of WIMP cold dark mattc;r searches 
is not good, since none of the many past or current ones 
has detected the cold dark matter particle ^But neither 
should they still do so if dark matter is just neutrino mat- 
ter. The benefit of neutrinos over other dark matter can- 
didates is that their existence is beyond any doubt. We 
have both questioned the reasons to rule them out and 
derived their mass. They can be dense, in the Abell 1689 
center a billion per 5 cc, but they are non-relativistic and 
annihilate each other too rarely to allow observation of the 
decay products in sky searches. 

* * * 

We thank Yu. Baryshev, C. Gibson, R. Schild, R. 
Balian, J. Smit and J. Kaastra for discussion. 
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